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Abstract: An analytical expression for Fermi-Dirac integrals of arbitrary 
order is presented, and its applicability in obtaining an analytical EOS of a 
degenerate non-relativistic Fermi gas is discussed to some extent. 

Introduction 

Degenerate Fermi gases occur in a variety of systems studied in astro- 
physics and "laboratory" physics. Examples of such systems include white 
dwarfs and neutron stars,planetary interiors, ordinary metals and organic con- 
ductors. 

The general form of the equation of state (EOS) of a non-relativistic 
degenerate Fermi gas has been determined more than half a century ago 
( for example, Chandrasekhar, 1939) as: 

4-7T 

n e = — (2m e k B Tf 2 F 1/2 ( V ) (1) 



All the symbols have their standardized meanings ,and -^1/2(7?) denotes a 
particular case of the Fermi-Dirac ( FD ) integrals of order n: 

*»W-J l + exp[/?(e-//)] [Z) 

where /(e) = e n ,n6R,/i denotes the chemical potential,/? is the inverse 
temperature, and r\ = /3/i. 
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The practical problem of the numerical evaluation of the FD integrals has 
been the object of several recent studies ( such as Cloutman,1989; Antia,1993; 
Miralles and Van Riper, 1995). The purpose of the present paper is to obtain 
an analytical approximation to the FD integrals of arbitrary order, and to 
apply it to the problem of the EOS of a degenerate Fermi gas. 

The calculations 

Taking ks = l,and introducing a change of variables by 

e-fi = Tz (3) 
eq.(2) can be transformed into the following form: 

F . m = T r = fm<u + T r n^-n, r, ) dz 

J-ti/T 1 + exp[z\ Jo Jo l+exp|zj 

(4) 

The difference of values of / in the last integral of eq.(4) can be developed 
into series as 

f(jM + Tz) - /(/, -Tz) = Y$- i-m K -^-f {n \») (5) 

n=0 n! 

Inserting eq.( 5 ) into eq.(4 ) and using the fact that 



/•OO J.OI — 1 

Jo 1 +exp[^] v ' v 7 v ' 



where r(a) and C( a ) denote the gamma function and Riemann's zeta func- 
tion,one gets the following final form of the FD integrals of arbitrary order: 

n=0 

This expression is a generalization to arbitrary order and number of terms 
of existing partial results ( for example, Landau and Lifchitz,1976 ). 



r n 00 f( n )(ij) 

F n (fo)= f(e)de+T^ J — p[l-(-l) B ]T n (l-2-»)r(n+l)C(n+l) (6) 
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Discussion 



Inserting n — 1/2 into eq.( 6 ) and limiting the sum to terms up to and 
including T 6 ,it follows that 

FyM = ^ 3/2 [l + ^(-) 2 + ^(-) 4 + |^(-) 2 ] (7) 
3 8/2 o4(J /i 3072 jj, 

This result is impractical,because the chemical potential itself is a function 

of T and the number density of the Fermi gas. It can be shown from eq.(l) 
that this function can be approximated by the following developement: 

2 rp 4 rp _6 T 7T 8 T 

The symbol Hq denotes the chemical potential of the electron gas at T = 

K,and it depends only on the particle number density.Inserting eq.( 8 ) into 
eq.(7) and developing into series in T,one gets 

3 48 /j, 5700 /io 14b fj,o 

Note that we have managed to express this FD integral as an explicite func- 
tion of the particle number density and the temperature. This result is a dis- 
tnct advantage over some previous numerical work ( such as Cloutman,1989 
), where the argument of this integral was left in the form /3/i,without tak- 
ing into account the number-density and temperature dependences of the 
chemical potential . 

How can eq.( 9 ) be applied to in obtaining an analytical EOS of a degen- 
erate Fermi gas? Inserting, at first, the known result for /i into eq.( 9 ),one 
would get an expression showing explicitely the dependence of a FD integral 
of the order 1/2 on the particle number density and temperature. Going a 
step further and inserting this result into eq.(l),one would get the required 
EOS, in the form of an equation relating the number densitytemperature and 
various known constants (such as h and ks ).This result has the form 

n e - K x Wn e T 312 - K 2 Wn; 5/3 T 11/2 - K 3 Wn~ 3 T 15/2 - K 4 Wn- 13/3 T 19/2 = 

(10) 
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where n e denotes the number density of the electrons and all the other sym- 
bols denote different combinations of known constants. The low temperature 
limit of this equation,obtained by developing it into series in T up to terms 
including T 9 , is solvable within S. Wolfram's MATHEMATICA 3 software 
package in a few minutes on a Pentium MMX/166 with 32 Mbytes RAM. 
The solution has the form 

n = — ' — << >> 

6V4 [[ i _ 3K.WT 3 / 2 + 3K 2 T 3 - K 3 W 3 T 9 / 2 

(11) 

where << .. >> denotes terms omitted due to space limitations. 

We have thus obtained an analytical form of the EOS of a degenerate 
Fermi gas. The volume per particle, needed in some applications, is simply the 
inverse density. Our result is approximate, in the sense that the developement 
in eq.( 9 ) is limited to a small number of terms. Increasing the number of 
terms and applying this EOS to real physical systems will be the subject of 
future work. 

Acknowledgement 

I am grateful to Dr.Gradimir Vujicic from the Vinca Institute of Nuclear 
Sciences for a discussion which motivated this work,and to Mr. Kristijan 
Lazic for the installation of the computing facility. 

References 

Antia,H.M.: Astrophys. J.Suppl.Ser.,84,101 ( 1993 ). 
Chandrasekhar,S.: Stellar Structure, Univ.of Chicago Press., Chicago 
( 1939 ). 

Cloutman,L.D.: Astrophys.J.Suppl.Ser.,71, 677 ( 1989 ). 

Landau, L.D. and Lifchitz,E.M.: Statistical Physics, Vol. l,Nauka, Moscow 

( 1976 ). 

Miralles,J.A.and Van Riper,K.A.: preprint |astro-ph/9509124j ( 1995 ). 



4 



